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Executive Summary 


Simpson Weather Associates, Inc. has developed LAWS 
Simulation Models (LSM). to evaluate the potential impact of 
global wind observations on the basic understanding of the 
earth’s atmosphere and on the predictive skills of current 
forecast models ( GCM and regional scale). Under previous 
contracts, SWA has developed two basic algorithms for use with 
simulated doppler lidar wind profilers. The first, a Shot 
Management Algorithm (SMA), controls timing and placement of 
lidar pulses. The second, a Multi-Paired Algorithm (MPA), 
extracts horizontal wind components from the unique lidar radial 
velocity observations. Fully integrated ’top to bottom 1 LAW'S 
Simulation Models for global and regional scale simulations were 
developed .during this contract period. The algorithm development 
incorporated the effects of aerosols, water vapor, clouds, 
terrain and atmospheric turbulence into the models. Other 
additions include a new satellite orbiter, signal processor, line 
of sight uncertainty model, new Multi-Paired Algorithm and wind 
error analysis code. An atmospheric wind field library containing 
control fields, meteorological fields, phenomena fields, and new 
European Center for Medium Range Weather Forecasting (ECMWF) data 
was also added. SWA has used the LSM to address some key LAWS 
issues and trades such as accuracy and interpretation of LAWS 
information, data density, signal strength, cloud obscuration and 
temporal data resolution. 

A synopsis of key work performed under this contract 
f o 1 lows : 

* LAWS Simulation Model Upgrades 

Designed and developed a new satellite orbiter that 
simulates any orbital inclination angle from pure 
equatorial to sun-synchronous polar. 

Developed a version of the Air Force Geophysical 
Laboratory ( AFGL ) LOWTRAN 7 model that provides water 
vapor attenuation and aerosol backscatter profiles on a 
1° X 1° lat/long grid. ECMWF data profiles are used to 
approximated the natural variability of the optical 
properties . 

Modeled shot scale atmospheric turbulence using a Von 
Karman turbulence technique. The model insures that the 
shot scale turbulence is consis tent with the inter-shot 
scale turbulence flow st rue tures that generate the 
shear and the finer scale turbulence. Turbulence due to 
the following phenomenas were added; convection, wind 
shear, mountain waves and jet streaks. 



Incorporated the baseline signal to noise equation into 
the LSM. 

Incorporated the baseline line of sight uncertainty 
equation into the LSM. 

Developed a Wind Field Generator Library that contains: 

Control fields such as divergence, vorticity, 
deformation and translation. 

Correlated meteorological fields from random generated 
fields . 

Phenomenas such as hurricanes, AVEVAS , mountain waves 
and jet streaks . 

ECMWF meteorological 1.875° X 1.875° lat/long 
prof iles . 

Developed a LAWS error analysis model that provides 
measurement errors, sample errors and 
representativeness. 

* Error Minimization Study 

SWA examined simulated spaced-based lidar wind errors 
due to line of sight measurement errors. The wind 
errors were computed for co-located and non co-located 
laser shots (10km and 70 km shot separation). 
Decreasing the shot separation by 60 km resulted in a 
12 to 20 % increase in the number of wind speed errors 
for the 0 to 1 m/s range. 

* Wind and Aerosol Inhomogeneities 

SWA examined sampling errors in the vicinity of wind 
and aerosol inhomogeneities. Two Situations were 
considered. First, the maritime boundary layer and 
second, an elevated temperature inversion within the 
troposphere. Simulated lidar measurement errors were 
found to be 5 to 10% due to coincident backscatter and 
wind speed gradients. These errors rival other errors 
expected f rom the current LAWS s ampl ing strategy and 
anticipated signal to noise. 

* Global Cloud Study 

Cirrus clouds will have a significant impact on the 
performance of LAWS as currently designed. SWA 
performed a literature search to obtain a 
reasonable estimate of the global distribution of 



cirrus clouds. We examined several summaries of 
satellite-based cloud climatologies and produced an 
expected LAWS performance chart as a function of 
latitudinal zones. The table was later upgraded to 
include the presence of low, middle and high clouds 
coupled with the availability of atmospheric aerosols, 
where available. 

* Signal to Noise sensitivity study 

SWA performed a sensitivity study on expected signal to 
noise (SNR) using the baseline SNR equation. The study 
computed SNR as a function of nadir scan angle and 
considered the effects of aerosol backscatter, 
molecular attenuation and satellite altitude. For the 
current LAWS orbital configuration (i.e., a satellite 
altitude of 705 km and a scan angle of 45 degrees) in a 
maritime atmosphere, the surface and midlevel SNR (db) 
was 12.23 and -0.81, respectively. 

* Observing System Simulation Experiments (OSSE) 

1) SWA is conducting an OSSE with Florida State 
University (FSU) using a very high resolution global 
spectral model ( Kr ishnamurti ) . SWA has provided 10 days 
of LAWS simulated wind data that along with World 
Weather Watch (WWW) data is being used as input winds 
to the FSU spectral model. Output from the assimilation 
is being compared to output using only WWW data as FSU 
wind input. 

2) SWA prepared for a regional scale OSSE using the 
LAMPS model on the CRAY-XMP at MSFC. Tim Miller, MSFC, 
is exploring assimilation of LAWS data into mesoscale 
models. This OSSE also serves the Shot Management 
Algorithm by providing a testbed for evaluating various 
scan patterns and shot densities. 

* Baseline Atmosphere for LAWS trade studies 

SWA produced a baseline atmosphere for the LAWS Science 
Team. The baseline atmosphere is a gridded data base 
with prescribed aerosol backscatter, molecular 
attenuation, wind and turbulence profiles in a 100X100 
km 3 volume. The aerosol backscatter profile is a 
composite based on ground-based lidar measurements 
taken at JPL and WPL . The attenuation profile is 
intended to be a severe mean maritime profile and is 
based on LOWTRAN 7 data. 

During the period of this contract, SWA personnel 
participated in 1 workshop, presented 2 papers, wrote 2 
additional papers and attended 4 LAWS oriented meetings. 
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Simulated space-based Doppler lidar 
performance in regions of backscatter inhomogeneities 

G.D. Emmitt and S.A. Wood 

Simpson Weather Associates, Inc. 
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ABSTRACT 

The prospect of obtaining directly measured winds on a global scale has raised questions 
about the expected quality of the lidar wind measurements and the potential 

to sampling^attems and line-of-sight impediments. Extensive cotn P u ^^ s ^ U J^ -n regions 
going to address these and other issues. One source of measurement bias is found in regio 
of the atmosphere where gradients in both lidar backscatter and the winds occur together. 

The potential biases that result are identified and their magnitudes estimated. 

1 . INTRODUCTION 

* — 

A space-based Doppler Lidar Atmospheric Wind Sounder (LAWS) has been proposed by NASA as 
a facility instrument for its Earth Observing System. 1 Hardware feasibility and data impact 
studies are on-going. 2 The uniaueness of a lidar wind measurement gives rise to many ques- 
tions regarding^ the^accuracy and interpretation of the information. Fundamental que: 
are related to the data density (limited primarily by laser lifetime and scan ra ) , 

Sdt^.1 resolution. A Lid.r simulation Kodul <£»>*“' a ^ t j£*JJS3 1 ,1 h . 
address some of these issues and to find ways to maximize the information content within 
current hardware configurations and performance constraints . 3,4 In Pigurel , the general 
density of LAWS data is shown with vectors representing data combined ln "ighest 

resolution mode (which is not necessarily the most accurate mode since no averaging is per 

formed) . 



Given our current expectations of global distributions of backscatter at 10 6 or 9.11pm, 
wf» anticipate useful LAWS data in the lowest 5 km over most of the globe an PP 

trite 03 The^LSM^is^generally^used^t^examine ' th^global^performance of various lidar s can- 

K J ? as 2 * 

sr jsasa-fS 1 ^ s 

interest has motivated an “order of magnitude study. 

We have chosen to look at two situations: (1) the marine boundary layer (below 500 

meters) and (2) elevated temperature inversions within the troposphere. “ ing rise 

to^larg^vertical gradiln^ln^ackscattefira layer where the typical wind profile also 



shows a strong change with height (Figure 2). In the vicinity of elevated temperature 
inversions , one often finds backscatter "spikes" (Figure 4) and wind velocity shears due to 
decoupling at the density interface. For both of these situations, the net effect is that 
when a weighted average of the winds within a lidar sample volume is obtained, errors are 
introduced in making height assignments of the velocity information. 

This weighted sampling is not unexpected and is common to other remote sensing systems. 
However, the magnitude of the errors is noteworthy as are some of the implications of the 
resulting biases to the computations of such quantities as heat and moisture fluxes. 

2. MARINE BOUNDARY LAYER 

The distribution of backscatter and winds in the marine boundary layer have been general- 
ized in Figure 2. The backscatter profiles result from a composite of the LOWTRAN 7 Navy 
Maritime Model as well as some special data sets compiled by surface based CO 2 lidars . The 
wind profile is the standard log (z/z 0 ) form and the surface roughness (z Q ) is taken to be 
0.1 meters (rough seas). 

In Figure 3, we have plotted the errors due to sampling volumes of different lengths. 

For a 500 meter pulse length the errors approach 10% at 250 meters. The dotted lines are 
the errors that would have occurred if there had been no gradients in Bio.g* While some of 
these differences could be corrected by accounting for such sampling related problems, the 
general bias is towards an underestimation of the wind shear and therefore the heat and 



Fiqure 2. Ideal representations of the vertical Figure 3. Comparison of lidar measure- 

distribution of backscatter (10.6 pm) ment errors (observed speeds- 

and wind speed in the lowest 2 km of actual speeds) for different 

the atmosphere above the ocean. Back- pulse lengths for o e 

scatter gradients are consistent with case of no backscatter gra- 

observations near the ocean’s surface. dients and the case shown m 

Figure 2. 

3. ELEVATED INVERSIONS 

Marine inversions, nocturnal inversions or cloud generated inversions can cause aerosol 
flux convergence and result in a high concentration of aerosols near the base of the 

temperature structure. Figure 4 shows how the winds respond to the inversion by accelerat- 

ing above it. Compared to the marine boundary layer case, the patterns of sampling errors 
are considerably different (Figure 5). Not only is the magnitude of the errors different 
but also the sense of the error. Without any backscatter structures the maximum lidar mea- 
surement error is an overestimate; with an assumed backscatter feature at the inversion, tne 
maximum errors are underestimates. It is noteworthy that the magnitudes of the extreme 
errors increase with pulse length. 






Figure 4. Ideal representations of the 

vertical distribution of back- 
scatter (10.6 ym) and wind speed 
in the vicinity of an elevated 
temperature inversions. 



Figure 5* Wind speed errors (observed 

wind speed-actual wind speed) 
for three different lidar 
pulse lengths applied to the 
profiles in Figure 4. 


4. SUMMARY AND CONCLUSIONS 

The maanitude of the lidar measurement errors (5% to 10%) due to coincident S and wind 
speed gradients rival other errors that are expected from the proposed LAWS sampling * tra_ 
tegy and anticipated Signal-to-Noise Ratio (SNR) . Advanced signal processing and 
computation algorithms should be able to reduce the magnitude of the errors , s J°" d er to 
3 and 5 However, the general biases towards lower wind speeds will 

correct'and must be addressed in the ongoing Observing System Simulation Experiments (0 
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i. Introduction 

A space-based Doppler Lidar Atmospheric Wind Sounder 
f LAWS) has been proposed by NASA as a facility instrument for its 
Eartt^ Observing 1 System (Hoi) (NASA, 

and data impact studies are on-going (Emmitt and Ho ^ st °"’ ; 

The uniqueness of a lidar wind measurement gives rise to many 

questions regarding the accuracy and intsl ' preta J' 0, J density 

Information. Fundamental questions are related tt th i data dens.ty 
(limited primarily by laser lifetime and scan rates), cloud 
obscuration andSempcral resolution. A Lidar 

fLSM) has been developed to address some of these issues ana v 
find ways to maximize the information content within i the * curre 
hardware configurations and performance constraints (Emmitt a 
Houston, 1986; Biibro and Emmitt, 1987). 

Given our current expectations of global distributions or 
backscatter at 10.6 or 9.11 pm, we anticipate userul LAWS data in 
the lowest 5 km over most of the globe and in upper regions of 
the troposphere where thin ciouds, volcanic dust, or ° the ^ ae ^, 
may concentrate. The LSM is generally used to ex ^' na . n ^ e D ?^^‘ s 
performance of various lidar scanning techniques sampimg P a « e 
and alaorithms In this paper we use the model to look at a suoser 
of circumstances where strong aerosol gn adients _occur ^jeg.ons f 
significant wind gradients, giving rise to measurement biases wh 
will require special interpretation. 

We have chosen to look at two situations: (1) the marine 

boundary layer and (2) elevated temperature inversions within the 
troposphere. In the marine layer one expects strong 9 r ^ ,ents 
airborne sea salt near the ocean surface giving r 'se to large 
vertical gradients in backscatter in a layer where the typical I win 
profile also shows a strong change with height ( F '9 U ^ ’ . , * 

the vicinity of elevated temperature inversions, one often finds 

backer ? spikes" (Figure 2) and wind velocity shears due » 

decoupling at the density interface. For both wTridfLlth in 

the net effect is that when a weighted average of the J'nds witn 
a lidar sample volume is obtained, errors are introduced in mak 9 
height assignments of the velocity information. 

This weighted sampling is not unexpected and J®*°^ m f h n e t0 
other remote sensing systems. However, the ma ^ d f the 

errors is noteworthy as are seme ot the ^ he at and 

resulting biases to the computations or such quantities as heat an 
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moisture fluxes. 


2. Marine Boundary Layer 



atmosphere above the ocean. 

(Bottom) Comparison of lidar measurement errors 

(observed speeds minus actual speeds) for 
different backscatter profiles including 
the case of no backscatter gradients. A 
500 meter pulse length was assumed. 

(zo) is taken to be 0.01 meters (rough seas). The composite 
backscatter profile results from several special data sets compiled 
by surface based C0 2 lidars. The Navy Maritime Profile was 
obtained from LOWTRAN 7 and the NASA Lidar data was taken 
from data supplied by NASA Langley Rsearch Center (Ed Browei ). 
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Figure 2. Ideal representations of the vertical 

distribution of backscatter (10.6 urn) and wind 
speed in the vicinity of an elevated temperature 
inversion. 

In Figure 1 (bottom), we have plotted the errors associated 
with sampling with a 500 meter pulse. The free stream velocity 
was chosen to be 10 m s _1 . The dotted lines are the errors that 
would have occurred if there had been no gradients in backscatter. 
The error of -3 m s _1 at 250 m results from simple linear 
averaging of a logrithmic profile from the surface to 500 meters. 

The solid lines are for errors compounded by backscatter profile 
weighting. The last two backscatter profiles produce .2 - .8 m s _1 
additional errors for the given wind profile. While some of these 
differences could be corrected by accounting for such sampling 
related problems, the general bias is towards an underestimation of 
the near surface wind speeds and therefore the heat and energy 
fluxes over the oceans. 


3. Elevated Inversions. 

Marine inversions, nocturnal inversions or cloud generated 
inversions can cause aerosol flux convergence and result in a high 
concentration of aerosols near the base of the temperature 
structure. Figure 2 shows schematically how the winds respond to 
the inversion by accelerating above it. Compared to the marine 
boundary layer case, the patterns of sampling errors are 
considerably different (Figure 3). Not only is the magnitude of the 
errors different but also the sense of the error. Without any 
backscatter structures the maximum lidar measurement error is an 
overestimate; with an assumed backscatter feature at the inversion, 
the maximum errors are underestimates. It is noteworthy that the 
magnitude of the extreme errors increases with pulse length. 
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4. Summary and Conclusions 

The magnitude of the lidar measurement errors (595 to 1095) 
due to coincident backscarter and wind speed gradients rival other 
errors that are expected from the proposed LAWS sampling strategy 
and anticipated Signal-to-Noise Ratio (SNR). Advanced 
signai processing and wind computation algorithms should be able to 
reduce the magnitude of the errors shown in Figures 1 and 3. 
However, the general biases towards lower wind speeds will be 
much harder to correct and must be addressed in the ongoing 



Figure 3. Wind speed errors (observed wind speed- 
actual wind speed) for three different 
lidar -pulse lengths applied to the profiles 
in Figure 2. 


5. References 

Bilbro, J.W. and G.D. Emmitt, 1987: Assessment of error sources 

for one component wind measurements with a space-based 
Doppler lidar, Optical Society of America’s Fourth 
Conference on Coherent Laser Radar: Technology and 

Applications, Aspen, CO, July. 

Emmitt, G.D. and S.H. Houston, 1986: Assessment of measurement 

error due to sampling perspective in the space-based 
Doppler lidar wind profiler, Second Conference on Satellite 
Meteorology/Remote Sensing and Applications, Williamsburg, 

VA, May. 

Emmitt, G.D. and S.H. Houston, 1987: Impact of a space-based 

Doppler lidar wind profiler on our knowledge of hurricanes 
and tropical meteorology, A.M.S. 17th Conference on 
Hurricanes and Tropical Meteorology, Miami, FL, April. 

NASA, 1987: LAWS Instrument Panel Report, Chairman, R.J. 

Curran, NASA Earth Observing System, Voi. Ilg, NASA 

Headquarters, Washington, D.C. ORIGINAL PAGE IS 

OF POOR QUALITY 




Advantages of Approximate 
Shot Coincidence with a Space-Based 
Doppler Lidar 

G.D. Emmitt 


Simpson Weather Associates, Inc. 
809 E. Jefferson Street 
Charlottesville, VA 22902 


A space— based Doppler lidar wind sounder (LAWS) is being 
proposed as an EOS facility instrument. On a polar orbiting 
platform, LAWS is expected to provide a global set of wind 
profiles throughout the troposphere. These profiles will be 
bounded below by the earth’s surface or opaque clouds and limited 
to those areas having sufficient aerosols and/or thin cirrus. 

LAWS measures a line-of-sight (LOS) component of the total 
average wind vector (u, v and w) within a cylindrical sample 
volume of radius 10—20 m and length 300—500 m. Laser lifetime, 
in addition to other engineering considerations, limits the 
density of samples to 1 per ~ 1200 km 2 . The widely spaced LOS 
components obtained from relatively small sample volumes must 
then be combined to estimate the horizontal wind speed. It is 
clear that much of the time some averaging will be necessary to 
obtain reliable wind estimates. Just how much averaging will 
depend upon (1) the SNR for the individual samples, (2) the 
variance in the wind field at scales on the order of the sample 
spacing, and (3) the users requirements. The poorer the signal 
strength for samples within a given area the more averaging will 
be needed to reduce measurement uncertainty. The greater the 
variance within the real wind field, the more averaging is 
required to get a representative measurement. The user may 
define the level of averaging by specifying a desired resolution 
volume (e.g. , 200 x 200 x 1 km 2 for GCM assimilation or 75 x 75 x 
.5 km 3 for mesoscale research). 

The baseline configuration for the LAWS employs a fixed scan 
angle ( " 45°), a fixed scan rate (~ 6 rpm) and a fixed pulse 
repetition frequency (10 Hz). The pattern of shots resulting 
from these baseline parameters has been shown in Emmitt (1985). 
Better management of these shots to extend laser lifetime and to 
optimize sampling distribution has also been explored. Shot 
management options include scanner/pulse scheduled programming to 
achieve near coincidence for forward and aft shot pairs. The 
advantages of shot pair coincidence are under study within the 
context of general LAWS science objectives and desired 
accuracies . 

The issue of shot coincidence is related primarily to the 
assessment of various algorithms that take a very limited number 
of radial (LOS) velocity measurements within a specified area and 
generally an estimate of the horizontal wind components. 
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Currently being considered are three basic ways that the 
observations can be combined to obtain estimates of the 
horizontal wind vector. The first and most obvious way is to 
define a resolution volume and then use a weighted least squares 
analysis with the LOS components and their direction cosines. 

The weighting coefficients could be the SNR for each LOS sample. 
The result is a single estimate of the average horizontal 
components for the resolution volume* 

A second method is to use a variational analysis scheme with 
a numerical weather model (e.g., GCM) . In this case the L 
components are assimilated into the model as radial wind 
measurements and the model parameters are adjusted to optimize 
the agreement with the lidar data (and other data) with the 
resulting model wind used as the best estimate of the actual wind 

vector. 


A third method is to combine the LOS measurements into pairs 
with each of the two shots having a different perspective on the 
wind flow. Each pair is used to compute the horizontal wind 
components. These pairs are then weighted depending upon SNR, 
location within the resolution volume and the shot geometry o 
the two shots in each pair. The weighted pairs are then averaged 
to obtain, not only an estimate of the volume average, but also a 
first order estimate of the wind structure within the resolution 
volume . 


In developing the third approach, which we call the Multi- 
Pair Algorithm (MPA), it was clear that there were advantages to 
having the two shots in a pair occur in close proximity (< 1 Km) 
to each other within a selected layer. When the SNR for each 
shot was high (> 5 dB) , one could obtain the highest resolution 
(* 50 km) product possible with LAWS. However, if necessary, the 
wind estimates for the pairs could be weighted and averaged to 
produce lower resolution data sets. 

While shot coincidence can be argued for the MPA approach 
from first principles (i.e., common volume sampling), it is not 
so obvious that the first two approaches mentioned above benefit 
from such shot management. This issue will be addressed in e 

presented paper along with the results of some simulations 
currently underway. The first and third methods are currently 
being evaluated to address the following questions . 


1) which method provides the best estimate of 
an area averaged wind profile? 

2) which method provides the most accurate wind 
profile for a subgrid scale location? 

3) what is the impact of managed shot-coincidence 
on each method? 
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A Reference Atmosphere for LAWS Trade Studies 
S. A, Wood and G. D. Emmitt 
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1. INTRODUCTION 

A space— based Doppler Li dar Atmospheric Wind Sounder (LAWS) has been 
proposed by NASA as a facility instrument for the NASA Earth Observing 
System. A LAWS Simulation Model (LSM) has been developed to assess the 
impact of a spaced-based Doppler lidar wind profiler on global and 
regional features. Hardware feasibility and data studies are on- 
going. 1 » 2,3 The uniqueness of global Lidar wind measurements from space 
r ai ses many fundamental questions that may impact the design of such a 
system. The distribution of aerosols that provide backscatter , the 
molecular attenuation that reduces signal strength, the effects of wind 
shear and turbulence that effect measurement accuracy, and the presence 
of thin cirrus clouds that can enhance the performance are all issues 
that must be considered. 

This paper describes a candidate reference atmosphere from the LSM*s 
atmospheric library. The reference atmosphere is used to e x a min e LAWS 
baseline signal-to-noise and line-of-sight velocity errors. 

2 . REFERENCE ATMOSPHERE 

The LSM atmospheric library provides a probabilistic aerosol 
backscatter profile, a probabilistic thin cirrus cloud backscatter 
profile, a molecular attenuation profile, a zig-zag wind shear profile, 
sub-pulse scale turbulence , and a correlated horizontal wind field within 
a 100 X 100 X 15 km 3 volume. 

The probabilistic aerosol backscatter profile, shown in Figure 1., 
was constructed from ground based lidar data taken at JPL and WPL. The 
circles indicate the median value (including data "dropouts" ) as a 
function of altitude. The number in the circles is the percentage of 
total observations associated with that particular median. The + 1 sigma 
error bars were computed from several hundred profiles. The model 
assumes that backscatter is log normal around the median at all levels. 

It is noted that the backscatter near the ocean surface is thought to be 
much higher than shown. The JPL and WPL lidar data sets did not have any 
contributions of thin cirrus clouds to the upper tropospheric 
backscatter. Therefore, the cirrus mode from 7 to 15 km has been 
estimated based on general reports high frequency of occurrence of 
thin subvisual cirrus clouds . The distribution of subvisual cirrus has 
been estimated as 30% at 7 km as seen from a ground perspective at JPL, 
Boulder and Hawaii and 50% at 14 km. This is believed to be 
underestimated for the tropics from a space perspective, where 70-80% may 
be the closer value. The cirrus relative backscatter is also assumed to 
be log normal. 
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The molecular attenuation profile, shown in Figure 2., was generated 
by LOWTRAN 7 model and represents attenuation in a tropical maritime 
atmosphere, Earth’s surface . No cirrus cloud attenuation is included. 

The atmospheric generator creates a "zig-zag” wind shear profile, as 
shown in Figure 3. This shear profile allows the effects of wind shear to 
be considered at any level in the atmosphere. A very general sub-pulse 
scale turbulence due to wind shear is included . Using Von Karman (-5/3) 
turbulence spectra for wind shear 4 , the LSM integrates the spectra over 
the pulse length scale, which is multiplied by an estimated total wind 
shear turbulence that is proportional to the "zig-zag” shear . 
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3 . A REFERENCE ATMOSPHERE APPLIED TO SNR AND LOS UNCERTAINTY 

The reference atmosphere ’ s median backs catter profile , a tropical, 
maritime attenuation profile and a shear layer of 0.005 s _1 was used to 
examine baseline LAWS signal to noise and line of sight velocity error , 
which is based upon pulse— pair autocorrelation processing of the Doppler 
signal. Figure 4. highlights that no SNR was near 5 db for the mid- 
levels nor at extreme scan angles at the surface. If 5 db is the 
threshold SNR for extracting useful line of sight wind measurements , then 
for a scan angle of 45 degrees, a backscatter greater than E-10 m-1 
sr-1 is needed. Figure 5. shows that the probability of getting the 
backscatter needed to obtain a 5 db SNR is nearly 80 % of the time at the 
surface, but quickly decreases to below 50% in the mid— level to 17% at 
upper levels. Sub-visual cirrus can increase the probability of getting 
5 db from 20 % around the tropopause to 50% at 14 km. 

If we could extract information at a lower threshold SNR, via some 
advance signal processing, the picture changes significantly. Figure 6. 
shows that the probability of getting backscatter to obtain a -5 db SNR 
is much higher, on the order of 80 % at the upper levels. Figure 7. shows 
the radial velocity uncertainty as a function of signal to noise. Errors 
on the order of 1 - 2 m/s are expected at the surface layer, where SNR is 
13 db. At a snr of 5 db, errors on the order of 8 m/s should be 
expected. Again, if an advance signal processing scheme could relax the 5 
db threshold by 10 db, then radial velocity errors at 5 db could be on 
the order of 1 m/s. 
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4. CONCLUSIONS 

We have defined one possible candidate reference atmosphere from the 
LAWS Simulation Model. We have looked at the baseline signal to noise 
and radial velocity error's using the reference atmosphere. Based on a 
tropical maritime atmosphere, we have shown that obtaining wind 
information in the mid-levels will be difficult unless better signal 
processing is possible and/or sub— visual cirrus is present. This study 
does not consider cloud obscuration, particularly in the PEL. A current 
follow on study is including clear line of sight cloud statistics for 
penetrating cloudy regions. 
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BASELINE ATMOSPHERE FOR PHASE I 
LAWS CONFIGURATION TRADE STUDIES 

Definition 

The baseline atmosphere is a gridded data base with 
prescribed backscatter, absorption, wind and turbulence profiles 
within a 100 x 100 x 20 km 3 volume. 

Purpose 

The. primary purpose for providing a baseline (or reference) 
atmosphere is to have a common basis for evaluating and comparing 
various LAWS lidar/scan configurations suggested by the two 
contractors. The intent is to generate a description of the 
atmosphere that will permit parametric expressions of LAWS system 
trades involving measurement accuracies (LOS and horizontal 
vector), representativeness, resolution and areal coverage. 

It is anticipated that the resulting expressions of critical 
cost and performance trades between resolution, coverage and 
accuracy will enable the LAWS Science Team to define higher order 
simulations required to select the optimum LAWS configuration to 
meet the science objectives. 

The baseline atmosphere is constructed with the following 
considerations and assumptions : 

1) Each contractor will have their own higher order 
atmospheric models for more detailed or focussed 
trade studies ; 

2) The LAWS Science Team needs to know the general 
sense of the costs ($ and performance) associated 
with various system configurations; 



3) The final LAWS configuration will be the 
product of several iterations between the 
engineering and science efforts; 

4) The most critical questions to be addressed 
by the baseline atmosphere are: 

* with what frequency will LAWS obtain 
accurate (+_ 1 m s - 1 ) winds along cloud- 
free line-of-sights (CFLOS) as a 
function of laser/optic parameter; scan 
angle; pulse length; pulse averaging; 
and height in the atmosphere* 

* what is the optimum set of baseline (+.) 
system parameters that will achieve the 
most accurate LOS measurements; the most 
representative LOS measurements ; the best 
global coverage; the most accurate horizontal 
wind estimates ( primat ive* ) ; the most repre- 
sentative horizontal wind estimates? 

^primitive - a Level 2b wind vector computed using only 
information provided by the LAWS instrument. 

Description 

The baseline atmosphere contains information on aerosol 
backscatter, absorption and winds applicable to a 10.6 jam or 9.11 
pm LAWS. Each component is treated independently while 
recognizing potential correlations in the real world. However, 



any attempt at this time to be more realistic is not justified. 
Therefore, the backscatter profiles, aerosol profiles, wind 
profiles and wind horizontal structures are designed to contain 
only information necessary to answer the critical questions 
listed above. 

Backscatter 

The GLOBE backscatter summary profile used in this baseline 
atmosphere is shown in Figure 1 . This format is attractive in 
that it highlights the mode of backscatter return with height and 
at the same time allows the likelihood of a specific backscatter 
value to be expressed in terms of the distribution. Both aerosol 
and thin cirrus profiles are presented. 

It is recognized that the GLOBE data set is based upon a 
limited body of observations and is probably a conservative 
representation of the backscatter over the entire globe. 

Absorption 

The water vapor and gas absorption profile was chosen to 
represent the tropical marine atmosphere. The arguments for this 
choice are that the current data sparse regions are the oceans 
(70% of the globe), the geostrophic approximation doesn’t apply 
well in the tropics (40% of the globe), and the absorption will 
be greatest in this region due to water vapor. Figure 2 is taken 


from LOWTRAN 7.. 



Winds 


In keeping with the generally schematic nature of the 
baseline atmosphere we have chosen to describe a 3-D wind field 
that will test the general performance of any algorithm to 
estimate the LOS and/or horizontal wind vectors. The primary 
attributes of the reference field are: 

1) vertical shear of the horizontal wind (selectable 
with 5 x 10" 3 s _1 as a default); 

2) horizontal coherent structure across the reference 
volume ( du/dx = dv/dy = 10~ 3 s - as a default) ; 

3) isotropic turbulence scales below 1 km and a 
vertical structure (ct 2 = 16 m 2 /s 2 to a 2 = 4 m 2 /s 2 ); 
correlated variance at grid scale (1 km); and 
vertical velocities not correlated from shot to 


4) 

5) 


shot - aw 2 


1 m/ s . 
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Subject: Contribution of shear to velocity spectrum 

broadening 

Given horizontal and vertical linear shears of the wind 
within a LAWS sample volume, what is the contribution to the 
spectral broadening? 

The following computations are based upon material presented 
in Doviak and Zrnic (1984), pp. 87-90. 

For LAWS (baseline), the following assumptions are made: 

x = pulse length in seconds = 4.0 ps 
CTd = Gaussian diameter of beam at 1000 km range = 10 m 
EL = scan elevation from horizontal = 37° 

K* = Ky = linear horizontal shear over 1 km = 10 " 3 s ‘ 1 
Kz = linear vertical shear over 20 km = 5 x 10 " 3 s - 1 
a 3 2 = variance in LAWS velocity spectrum due to linear 
shear across the sample volume 

Os 2 = ( CTd Ke ) 2 + (od K 0 ) 2 + (or Kr ) 2 

ar = .35 c t/ 2 = 210 m 

0,0 = angular direction around LOS, r0 is parallel to 
the ground 

An order of magnitude argument yields: 
as 2 = (10 x 10- 3 ) 2 + (10 x lO " 3 ) 2 + (210 x 10- 3 ) 3 
= 10 ' 4 + 10 * 4 + 4.4 x 10 “ 2 

as 2 • (Or Kr 2 ) 

Kr = (Kz 2 + Ky 2 ) 1 / 2 cos -EL + Kz • sin EL 
= 1.12 x lO " 3 + 3.0 x 10 * 3 = 4.12 x 10 ' 3 
da 2 = (210 * 4.12 x 10 - 3 ) 2 = .748 m 2 /s 2 


Thus : 
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LAWS OSSE Using FSU 

Global Spectral Model 
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An illustration of the transform grid at resolutions 
T21, T1Q6 and T170 over a part of the globe. (Triangular 
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SWA LAWS Simulation of 
Laser Pulses Intercepting 
Clouds in the FSU OSSE 


* Use ISCCP Data 

— Optical Depths < 1.3 are Cirrus 
— Atmospheric Layers (0— lOOOmb) 

— % Coverage (2.5 X 2.5 deg) 

* At Each Level Calculate Probability of 

- Clear (No Clouds) 

— Thin Cirrus 

— Opaque Cloud 

* Use ISCCP Data to Provide Probability 
of LAIS Receiving Backscatter from 
Sub-vis ibile Cirrus 













Preparation for OSSE Using a 
Regional Scale Model — LAMPS 


Tim Miller 

NASA Marshall Space Flight Center 
G.D. Emmitt 

Simpson Weather Associates 


Objectives : 

- Develop strategies for assimulating LAWS data into 
a non-hydros tatic model - line-of-sight , horizontal 
vectors, hybrids, etc. 

- Assess impact of LAWS data on model forecasts 

- Provide testbed for evaluating various LAWS scan 
patterns and shot densities as part of the 

Shot Management Algorithm. 
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A. Monthly Act ivi ties 

During the period 9/ 1 9*“ 1 0/ 1 9 we returned to a full effort 
following a 6 month period without new funding. In the last 
month we have 

1) Incorporated a modified SNR into our LAWS simulation 
model ; 

2) Using LOWTRAN 6 and the LAWS SNR we modelled the effects 
of subvisible cirrus on the velocity estimation error; 

3) Began development of a tri-scale wind field model to be 
used in the LAWS Analysis Reference Wind Fields; and 

4) Defined a new data quality function for use in the 
Goddard GCM OSSEs, 

B. Next Month’s Work Plan 

We will be focusing upon the tri-scale wind model with the 
goal of having a preliminary set of reference fields ready by the 
end of December. We will also prepare a bi-monthly technical 
report with details on our simulation activities. 
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A. Current Status 

During the period 10/19/88-11/18/88 we completed the installa- 
tion and checkout of LOWTRAN 7 absorption and backscatter profiles for 
the entire globe on a l°xl° lat/long grid. These profiles will be. 
used in the computation of the SNR of individual lidar shots used in 
conducting the OSSEs. 

We are continuing the development of the "tri-scale" wind variance 
model for use in the LAWS Analysis Reference Wind Fields. Our basic 
approach is to use the Von Karman model for several wind phenomena such 
as mountain waves, jet streaks, tropical storms and dry convective 
boundary layers. 

A technical report on the LAWS effort since June 1988 is being 
sent under separate cover since it reports on related work not funded 
under this contract. 


B. Over the next two months we will continue to evaluate the "tri-scale" 
model as well as improving the LSM for execution on the PC-AT. We 
are also preparing a paper for presentation at the January meeting 
of SPIE. 
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A. During the last month, progress has been made in developing 
a fully programmable wind variance model covering the range 
of a few meters to several 100 kilometers. The purpose of 
is to be able to run Monte Carlo tests for LAWS sampling of 
various mesoscale phenomena such as jet streaks, mountain 
waves, tropical storms, etc. The variance model insures 
that wind variance on the lidar shot scale is consistent 
with the inter-shot scale flow structures generating the 
shear and thus the finer scdie turbulence. 

A paper has been written and submitted to SPIE for inclusion 
in the proceedings of the SPIE’s OE/LASE ’89 in January 1989. 
A copy is included in Appendix A. 

B. During January, we hope to complete the prototype for the 
LAWS Reference Wind Field Package to be used in conducting 
system trades during the upcoming feasibility studies. 
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A. Current Status 


We are continuing the development and Implementation of the TRI 
scale" wind variance model into the LAWS simulation model. It is 
recognized that the wind variance on a given scale depends on the state 
of organization of the atmosphere and the use of a single value to 
represent variance is unrealistic. Our current approach is to use the 
Von Karman model to provide spatial and temporal estimates of natural 
wind variance for our LAWS analysis reference wind fields. 


We have selected to model four basic turbulence fields: convec- 

tive processes, wind shear, rotor waves and mountain waves based on 
experimental data described in Rhyne et al. (1976). Rhyne used tur u- 
lence sampling programs to describe each turbulence field using a power 
spectral density function, the Von Karman model (equation 1). 


<J> (1/X) = 2a 2 L 


fl + 8/30.339 L 2 ttA) 2 ] 
[1 + (1 .339L 2ttA) 2 ] ,,/6 


( 1 ) 


In this equation <5 is the power, a 2 is variance, A is wavelength and L 
is the length scale of the turbulent phenomenon. Figure 1 shows a curve 
generated by the Von Karman equation depicting the variance at long 
wavelengths and at shorter wavelengths within the inertial subrange for 
different scales of organized turbulence. 


We partitioned the Von Karman wind variance into three parts 
defined by the LAWS sampling geometry. The three divisions are based on 
three scales of sampling: (1) Shot Scale (SS) , 1000 meters or ess. 

Defined from the major axis of the cylindrical pulse volume; (2) Pulse 
Scale (PS), 1000 to 100,000 meters. Defined based on spacing between 
shot pairs for the estimation of wind; and (3) Large Scale (LS) , greater 
than 100,000 meters, based on the stated goal of LAWS to provide wind 
profiles with a horizontal resolution of 100 km. 

Using these scales, the area under the power density curve can be 
divided into three sections with the "knee" of the curve fa 1 1 i ng wi th . n 
the PS area. When determining the variance, we need only determine the 
area for SS and LS and by subtraction from one (the total normalized 
variance), we can determine the PS value of variance. We reduced equation 
(1) in the following manner to determine the percent variance in each ot 
the three sampling ranges: 


${1 A) = 2a 2 L 8/3 (1.339 L 2 ttA)" 5/3 

= (2 * 8/3 * .0287) * a 2 * L _2/3 * (1A) 5/3 
= .153 a 2 L" 2/3 (1A)“ 5/3 


The area in the SS region is: 



2 


A (SS) = 


CO 

J <i> (1/X) * dd/X) 


.23 a 2 L - 273 * ( 1 A )" 2/3 



In the LR area Von Karman's equation reduces to: 
$(1/X) = 2 a 2 L 


I integrating : 


A ^LR^ 


IQ' 5 


$(1/X) d ( 1 /X ) 


>10 


= 2o Z L J d(l/X) 
o 

= (2*L*10~ 5 ) * a 2 

2 

The a 2 PS area we get by subtraction: a PS = 


, 2 2 
1 - a LR + ° SS 


The distribution of the total variance between the three areas can 
>e substantially different for different atmospheric phenomena With 
set bounds on the scale of the three divisions, .t is the length scale of 
the turbulent phenomenon itself, L, which determines the apport .onment 
of variance under the curve. 

Rhyne et al. (1976) measured all three components (longitudinal , 
lateral and vertical) of the wind for the four classes of J u, ^" C ®; ales 
Their data fit to the Von Karman power spectra at different 9 
and their values for standard deviation for the three components ^uv 
and w) and the appropriate length scales are given in Table 1. To f. 
the percent of total variance, we integrated equation (1) and based o 
the experimental standard deviation values, computed the actual variance, 
a 2 , for the four classes of turbulence. 

In a general sense, using the proven relationship be 5? ee " ^ . 

turbulence phenomena and the Von Karman model , we are now able to simulate 
total variance at three scales (l.e., shot scale, pulse scale and large 
scale). Next month we will be running LSM situations on °ur LAWS anal ys 
reference wind fields with the Von Karman turbulence rep ^ se "!^^ Variance 
are also continuing our development of a fully programmable . 

model that will allow us to simulate LAWS sampling of various mesoscale 

phenomena . 


Reference: 


Rhyne, R.H. 

^er Spectral 

Meteorological Cond I tions 11 


and K. Sidwel 1 , 1976: 


H.N. Murrow, aru rv. ■ ----- • , 

Turbulence Power Spectral Measurements to Long Wavelengths ^Severa 


NASA 


Conference, Hampton, VA, pp 


'Atmospher i c 
Severa 
271-286. 



3 


B. A paper was presented at the SPIE's OE/LASE '89 Anaheim, CA 
conference . 









Table 1. Four cases of turbulent phenomena from Rhyne et al. (1976). 
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George D. Emmitt 
S . A . Wood 
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A. During the period January 19-February 18, 1989, we have 

continued our development of the M LAWS Reference Fields” 
and the investigation of LAWS data interpretation in the 
vicinity of pulse scale inhomogeneities in the winds and 
backscatter , In particular, we have been simulating the 
advantages of achieving shot coincidence ( intersecting 
pairs) compared to the improvements obtained with higher 
SNR for LOS measurements (see attached report). 

In January, a paper entitled "Simulated Space-Based Doppler 
Lidar Performance in Regions of Backscatter Inhomogeneities" 
was presented by G.D. Emmitt at the 1988 SPIE conference in 
Los Angeles. A follow -on paper by Emmitt and Wood has been 
prepared and submitted for presentation at the 5th Conference 
on Coherent Laser Rdars to be held during June 1988 in 
Munich, FRG (see Appendix A). 

B. Most of our effort for the next month will be focussed upon 
the documentat ion of the LSM, simulation of the current set 
of reference wind fields, and further investigation of a 
hybrid approach to the use of LAWS data, i.e., use of the 
LOS measurements by themselves (regressed against a model 
wind field) and/or shot pair estimates when available. 



Error 


Minimization Study 


The Lidar Simulation Model ( LSM ) was run for a pure 
convergent field with and without the Von Karman modeled 
convective turbulence. Figure 1 shows the error field produced 
from the input wind field and the lidar simulated wind field with 
random perturbations in the line of sight measurement (1 m/s). 

The shots were assumed to be co-located. Wind errors greater 
than 1 m/s have been highlighted in red (in original document). 
Figure 2 is the same as Figure 1 except that only random 
perturbations due to convective turbulence have been considered. 
In this second case the lidar shots were not co-located and had a 
separation distance of approximately 70 km. Both figures show a 
significant number of wind errors greater than 1 m/s. 

The distribution of the wind errors was further examined for 
various scales of line of sight errors and convective turbulence, 
i.e., shot separation. Figure 3 is a histogram of the number of 
wind errors (u and v components and wind speed) falling within 1 
m/s increment error bins. The error bins go from 0 to 15 m/s and 
the last bin represents errors 15 m/s and greater. The top graph 
is for a line of sight error of 1 m/s and the bottom graph is for 
a line of sight error of 0.5 m/s. The figure shows a shift to 
lower (47%) wind speed errors by decreasing the line of sight 
error . 


Figures 4 and 5 are histograms of wind s 
errors for a convective turbulence case with 
10 and 70 km. The top graph as a LOS error o 
graph has a LOS error of 0.5 m/ s and the bott 
error. Table 1 summarizes the approximate pe 
the number of wind speed errors in bin 0-1 m/ 
the shot separation by 60 km. The percentage 
12-20% for the three LOS errors. 


peed measurement 
shot separations at 
f 1 m/s, the middle 
om graph has no LOS 
rcentage increase of 
s due to decreasing 
increase is around 


Figures 3-5 and Table 
has the potential of signi 
errors, the sampling relat 
general result lends itsel 
SNR improvements vs active 
coincidence . 


1 suggest that while improving the SNR 
ficantly decreasing the horizontal wind 
ed errors mask that advantage. This 
f to a cos t/bene? f i t analysis involving 
shot management (required for 



List of Figures 


Figure 1: Simulated space-based lidar wind errors with 

random perturbations in the line of sight 
measurement (1 m/ s ) . The shots in shot-pair 
are co-located. 

Figure 2: Simulated space-based lidar wind errors with 

random pertubations due to convective 
turbulence. The lidar shots are not co- 
located. No line of sight measurement error 
was considered. 

Figure 3: A number distribution of simulated space- 

based 1 i dar wind errors for tw r o line of sight 
errors. No atmospheric convective turbulence 
was considered. The shots are co-located. 

Figure 4: A number distribution of simulated space-based 

lidar wind errors for three line of sight errors. 
Atmospheric convective turbulence was considered. 
The shots are approximately 10 km apart. 

Figure 5: A number distribution of simulated space-based 

lidar wind errors for three line of sight errors. 
Atmospheric convective turbulence was considered. 
The shots are approximately 70 km apart. 
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A. During the period 18 February to 19 March, we 


1) continued to construct and evaluate reference wind 
fields for use in hardware trades and LAWS feasibility 
demonstrations . 

2) updated cirrus global climatology with recently 
published reports * 

3) prepared for and attended the EOS meeting in 
Washington, D.C. 

4) prepared for a "bi-modal" OSSE using ECMWF nature 
runs , and 

5) began preparation of a ’'LAWS scan geometry” packet 
for members of the LAWS facility team. 

B. From 19 March to present we have been working on a no-cost 
extension. Anticipating a new contrct being in place by 7 
April, we will begin focusing most of our effort on conduct- 
ing the Global OSSE and preparing for some regional scale 
model OSSEs with Dr. Miller (NASA/MSFC) and with Dr. 
Krishnamurti (FSU) . 



Expansion on A. 2 


Cirrus clouds could have a potentially significant impact on 
the performance of LAWS as currently designed. Dr. M. Hardesty 
(NOAA, WPL ) has proposed (as a member of the LAWS team) to 
examine the performance of LAWS in regions of cirrus clouds. In 
our case, we are interested in obtaining reasonable estimates of 
the global distribution of cirrus as input to our efforts to 
perform OSSEs ♦ 

A review of the most current published data and papers on 
global cloud distributions has revealed both useful yet ambiguous 
data - useful, in that the general distribution of clouds by 
season are reasonably represented using combinations of ground- 
based and satellite observations; ambiguous, in that the 
distinction between thin cirrus and opaque cirrus is vague enough 
to support a broad range of assumptions. 

The accompanying figures represent the summaries of several 
satellite-based cloud climatologies. There are several notable 
features as far as LAWS is concerned: 


1) A generally tri-modal distribution of cloud cover 
fraction between 70S and TON with peaks near 60-80% 
at 60S, 5N and 60N (minimums of 30-40% at 20S and 30N) 
[ Figures 1 and 2 ] . 


2 ) 


Conflicting 
from 40-90% 
1 ] . 


reports of cloud cover 
(seasonality may be an 


near the poles ranging 
explanation) [Figure 


3) High 10-14 km) clouds (opaque?) show a similar 
distribution with latitude - i.e., tri-modal but a 
lower fractional coverage (^ 10-30%) [Figure 3]. 


4) Cirrus (transparent?) coverage ranges from 30-75% 
with the maxima at 50S, 5N and SON [Figure 4]. 


Figure 5 is a very preliminary attempt to integrate the 
information in Figures 1-4. The assumptions are that the high 
clouds in Barton’s (1983) study (Figure 3) are opaque to 10.6 m 
and that both the cirrus and thin cirrus in the SAGE study 
(Woodbury and McCormick, 1986) are transparent. 

The interpretation of Figure 5 (TROPICS) is as follows: 


1) 90-100% of the time LAWS will get a return from the upper 
troposphere because of the presence of either opaque cloud or 
transparent cirrus. The return from aerosols at those levels is 
left as unknown at this time but expected to be certainly non- 
zero . 



2) 40-50% of the time there will be clouds in the mid- 
troposphere from which LAWS will get a return. The degree to 
which aerosols will provide sufficient backscatter during the 
remaining 50-60% of the time is an unknown , 

3) By definition, there are no clouds in the PBL and thus 
LAWS will get returns there only from aerosols. However, in the 
tropics we can expect to see into the PBL only 20-30% of the 
time . 

While Figure 5 presents the general picture for LAWS 
returns, there are seasonal variations in cloud cover that would 
need to be incorporated into any comprehensive impact assessment. 
In the short term, these results can be used to address seme 
basic questions regarding the "value" of PBL winds vs upper 
tropospheric winds in GCMs . We will continue to update and 
refine the interpretation of the available cloud climatologies, 
in particular those derived from SAGE II observations. 

It should be noted that any estimate of thin cirrus coverage 
is likely to be an underestimate due to detection thresholds by 
current space- based 


sensors . 



TOTAL CLOUD FRACTION (%) 



LATITUDE 


Figure 1. Zonal ly averaged total monthly mean cloud amount for July. 

The Nimbus (N7GCC) TH I R infrared noon and midnight cloud fraction 
estimates are compared with the N7GCC bispectral noon-cloud- 
fraction estimates for July 1973- Also shown in the figure are 
results from Beryland and Strokina's (1980) 30-year cloud 
climatology (BSS) , the Air Force three-d imens iona 1 -nephana lys i s 
cloud data for 1979 (AF 3D-N) ; compiled by Hughes and Henderson- 
Sellers, 1985; and London's (1957) multiyear averaged (northern 
hemisphere) cloud climatologies. Mote that London's southern 
hemisphere cloud amount is taken from the value in the northern 
hemisphere for the opposite season. (After Hwang, P.H. et al . , 
1988: The Nimbus-7 Global Cloud Climatology, BAMS , 7^3-752. ) 





Figure 2. Zonal ly averaged 3D-nephana 1 ys i s cloud amount for January, 
April, July and October 1379 * (After Hughes, N.A. and A. 
Henderson-Se 1 lers, 1385: Global 3D-Nephanalys is of Total 

Cloud Amount: Climatology for 1979 * JCAM , 669 * 686 . 




LAT1TU0E 


Figure 3. Seasonal zonal average distributions of the occurrence of 
high clouds with D1 reflectance greater than 0.01. (After 
Barton, I.J., 1983: Upper level cloud climatology from an 

orbiting satellite, JAS, •) 
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Expansion on A. 5 


The LAWS facility team is composed of individuals familiar 
with past LAWS simulations as well as new people who have asked 
for a brief summary of LAWS scan geometry and spatial coverage. 
The following material will be included in a status paper on LAW T S 
shot management being prepared for the LAWS team* 


The perspective of LAWS on the earth’s atmosphere is 
sketched in Figure 6* The values for the pertinent angles and 
distances are presented in Table 1 for several space platform 
altitudes . 


PHI is the nadir scan angle at the satellite. 


THETA is the angle to the horizon of the lidar beam at the 
earth’s surface. Theta is always less than the compliment of the 
nadir scan angle PHI . 


COS represents the percentage of the horizontal wind 
component sensed along the line-of-sight (LOS). 

SIN represents the percentage of the vertical wind 
component sensed along the line-of-sight (LOS). 

SWATH is the width (km) of the scan pattern on the ground. 
For the polar platforms, the spacing between satellite passes 
will be ^ 2800 km at the equator. 

RANGE is the LOS distance (km) to the earth’s surface. 

SNR is the signal- to-noise ratio near the earth’s surface 
using a backscatter value of 3 x 10" 7 m" 1 str" 1 . 


The trades between nadir angle, global coverage (SWATH) and 
SNR are clearly challenging. Better coverage and better 
sensitivity to the horizontal speeds are tied to available SNR 
and, in a sense, oppose the desire for greater shot density 
(given a fixed PRF ) . 
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A. During the past reporting period we have continued to work on 
the following LAWS issues: 

1) Review of recent developments in Signal- to— Noise Ratio 
(SNR) equations and the Line-of -Sight (LOS) velocity 
uncertainty estimators ; 

2) Vertical distributions of atmospheric turbulence for 
10 m < A < 10,000 m; 

3) Cirrus climatologies derived from several satellite 
based cloud climatologies; 

4) Preparation of a LAWS team "white paper" on the 
subject of shot management; 

5) Preparation for the Munich conference on Coherent 
Doppler lidars; and 

6) Attendance at the GLOBE meeting in Huntsville, AL . 

A review and computational comparison of SNR equations and 
LOS velocity estimates has raised several questions: 


(1) What is the appropriate SNR equation for single shot 
lidar wind measurement? 

(2) What is the proper estimation of uncertainty in the 
LOS velocity? 


and 

(3) 


If the performance of the LAWS is being overestimated 
by 5-10 db ( Kavaya - personal communication) then 
there is cause for some concern regarding making 
measurements from a POP in regions where 3 < 10“ 9 
m* 1 sr“ x . In Table 1, the following equation was 
used to compute the SNR for LAWS for various 
platform and scan angle combinations. 


SNR = — 

n i 

where 


7T 

i s 

n 

is 

T 

is 

D 

is 

c 

is 

J 

is 

K 

is 

hv 

is 


J c T 


D 2 K 0 


-2/a dl/slmj) 


8 hv R (NF) 


a constant [3.14159] 
the system efficiency [0.1] 
the pulse duration [6.77 E 
the telescope diameter [1.5] 
the speed of light [3.0 E 8] 
the power transmitted [10] (J) 


( s ) 


- 6 ] 

( m) 

( m/s ) 


( J) 


a is the molecular attenuation coefficient (m _1 ) 

1 is the path length (m) 

$ is the scan angle (radians) 

R is the range (m) 

(3 is the backscatter coefficient (m _1 sr‘ 1 ) 

NF is hard coded noise factor due to speckle, 
j i tter , etc . [1.5] 

The noise factor is believed to be underestimated. 

From Table 1, it can be seen that for a 705 km polar 
orbit and a scan angle (PHI) of 45° that: 

- the viewing angling at the earth’s surface 
will be 38.25° from the horizontal. 

- 79% of the u and v wind components will be 
projected into the LOS. 

- 62% of the vertical wind component will 
be projected upon the LOS. 

- the width of a data swath will be 1501 km. 

• the LOS range from the lidar to the earth’s 
surface will be 1498 km. 

- the SNR for the boundary layer with a (3 of 
10" 7 m* 1 sr“ 1 will be 21.96 dB . 

- the SNR for the mid- tropospheric layer with 
a (3 of lO -10 m“ 1 sr" 1 will be 1.18. 

- there will be a 70% coverage of the globe’s 
surface in 12 hours. 

- there will be an 18% overlap of samples taken 
in 12 hours (mostly near the poles), and 

- there will be an average of 7 shots into a 
100 x 100 km area. 

Figure 1, illustrates the relationship between scan angle 
angle and SNR in the tropical PBL as well as the percent 
global coverage during a 12 hour period. 

If 5 dB is the minimum SNR for extracting a useful estimate 
of the LOS wind component and the equations overestimate 
the M real " SNR by 5 dB, then good measurements in the 
mid-troposphere will require a § greater than 10‘ 9 m* 1 
sr - 1 . 


B. Effort will continue on the above items over the next month. 
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A. During June our efforts have been focused upon the 

development of a reference atmosphere for use by the Phase 
A/B contractors ( GE and Lockheed) in their initial trade 
studies. The motivations for having an '’initial trades” 
reference atmosphere are: 

1) to provide a reasonable "link" between line-of-sight 
measurements and the estimates of horizontal wind 
components in a turbulent and inhomogeneous atmosphere; 

2) to allow preliminary hardware configuration tracks 
to be explored prior to a full evaluation using more 
complete and globally representative simulations; 

3) to allow comparisons to be made between GE and 
Lockheed trade studies. 

Following discussions with M.J. Post (NOAA) and D. Bowdle 
(USRA/MSFC), we have decided to generate a "single situation" 
model that includes the most critical environmental elements 
for simulating a LAWS measurement - except, that is, for 
clouds. Some of the questions that can be addressed with 
the reference field are: 

1) What is the contribution of pulse scale turbulence 
to the LOS wind estimate? 

2) How critical is the choice of pulse length with regards 
tp SNR, wind shear (non-linear), aerosol inhomogeneities, 
e t c . ? 

3) What is the most likely vertical distribution of 
velocity accuracy? 

4) How does the LOS velocity estimate change with 
scan angle? 

5) How representative are the LOS samples? 

6) What order of accuracy can be expected for the horizontal 
wind estimate? 

7) What are the trades between LOS accuracies and 
sample density? 

The reference fields are not designed to answer the broader 
questions of global performance with cloud contamination, 
optimal configuration for GCM impacts, etc. 

The goal is to deliver the reference fields by 7 Julv 
1989 . 



Two Observing System Simulation Experiments (OSSE) have 
been initiated within the last month. One is with Florida 
State University ( Kr ishnamurt i ) and the other with Marshall 
Space Flight Center (Tim Miller). 

The OSSE at FSU is being conducted using a very high 
resolution global spectral model developed by Kr ishnamurt i . 

S. Houston (SWA), is working on-site to install the LAWS 
Simulation Model on the MSFC’s EADS and to obtain a 
simulated LAWS wind data set from an ECMWF analyses. A 
10-day forecast impact will be evaluated with the spectral 
model. The simulation will include: 

clouds - derived from satellite observations for case 
study 

aerosols - profiles taken from LOWTRAN and modified with 
GLOBE findings 

water vapor - output from model 
winds - grid scale from model 

turbulence - parameterized from grid scale winds and 
gradients . 

The second OSSE is being conducted using the regional scale 
model, LAMPS, on the CRAY-XMP at MSFC. Tim Miller is 
beginning to explore issues regarding the assimilation of 
LAWS data into mesoscale models. A 3-D gridded field (LAMPS 
output) has been received at SWA. Simulations will begin 
within the month of July. 

During June, Emmitt attended the LAWS feasibility study 
Requirements Review held in Huntsville, AL . One key issue 
raised during those meetings with GE and Lockheed was that 
of the reference atmosphere coupled with a baseline 
configuration. As mentioned earlier, SWA is preparing a 
"quick look" reference atmosphere for preliminary trades. 

A presentation was made at the 5th CLRC in Munich, FRG - 
"Simulation of a Space-Based Doppler Lidar Wind Sounder - 
Sampling Errors in the Vicinity of Wind and Aerosol 
Inhomogeneities" (Emmitt and Wood). An additional 
presentation was made at the special session on SNR 
equations . 

During the month of July we plan to deliver a reference 
atmosphere to MSFC for use by GE and Lockheed and to continue 
work on the two OSSEs. 
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A. 


We continue to work on the baseline {or reference) 
atmosphere for use by the Phase I/II contractors. Attempts to 
incorporate GLOBE backscatter data have raised some issues that 
need to be resolved before delivering the reference atmosphere. 

Work also continues on preparing an OSSE at Florida State 
University. The global scale OSSE, proposed by Dr. Krishnamurti 
at Florida State University (FSU) to simulate LAWS, will include 
Simpson Weather Associates 1 (SWA) LAW 7 S Simulation Model (LSM) 
output winds plus the World Weather Watch (WWW) data as input to 
the FSU Global Spectral Model for a 10-day simulation beginning 
with 12Z July 5, 1984. The output from this forecast will be 
compared with the forecast derived from using only the WWW data 
as input to the same Global Spectral Model. It was decided that 
SWA’s polar orbiting LSM would use ECMWF w r ind data as input every 
12 hours beginning 12Z, July 5, 1983. The ECMWF u and v wind 

component data are gridded on 1.875° latitude x 1.875° longitude 
and are available at the mandatory pressure levels. The data set 
also includes geopotential height, relative humidity and 
temperature at each of the mandatory levels. Topography data 
will be provided by FSU for the model to determine the lowest 
possible surface layer of wind output in the absence of clouds. 
Because the LSM requires the cloud inputs to simulate the lack of 
laser penetration in disturbed moist regions of the atmosphere, 
FSU will provide cloud heights and optical depths every 12 hours 
from the International Satellite Cloud Climatology Project 
(ISCCP). FSU will also provide input support for the 
parameterization of atmospheric turbulence in the model. Aerosol 
molecular optical properties will be estimated by the LSM to 
provide backscatter and attenuation effects on the simulated 
laser pulses. Relative humidity data from the ECMWF fields will 
be among the parameters to be used to estimate the aerosol 
distributions and natural variability. Sample output from the 
LSM for 12Z, July 5, 1983 will be provided to FSU to verify that 

the ECMWF data were input correctly. 

A version of the LSM model has been transferred to the 
Engineering Analysis Data System (EADS) at Marshall Space Flight 
Center (MSFC). The model is now running with input from the 
ECMWF u and v wind component, geopotential height data and 
gridded topography data for the globe. Work is also underway to 
develop a means to include atmospheric turbulence in the LSM. 

LAWS was represented by G.D. Emmitt at the recent EOSDIS 
Architecture Review held GSFC. A report on that review and 
actions taken by the EOSDIS Science Advisory Panel is being 
prepared for presentation at the LAWS Science Team meeting 
scheduled for August 10-11. 

B. 


We will continue work on the OSSE and baseline atmosphere. 
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A. 


During the last month we have made progress on several 
tasks : 

1) incorporated latest backscatter profiles into the 
LSM; 

2) developed a Baseline Atmosphere for initial trade 
studies on LAWS configurations; (see Appendix 

A) ; 

3) developed (but not tested) an algorithm for 
propagating LAW’S shots through a cloudy 
atmosphere using ISCCP data; (for use in OSSE’s); 

4) attended and presented material at the 2nd LAWS 
Science Panel Meeting (9-11 August 1989) at 
Huntsville, AL ; 

5) attended (as an EosDIS Science Advisory Panel 
Member) the EosDIS Architecture Review meeting 
at GSFC (24-28 July 1989). 

B. 

In the remaining 2 months we intend to complete the Baseline 
Atmosphere checkout and provide a statistical module for data 
retrievals based upon Baseline (3 profile (with and without 
cirrus ) . 
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A. 


During this reporting period we have: 

1) Delivered a LAWS Baseline Backscatter and Absorbtion 
profile and a baseline SNR equation for use by the 
Phase A/B contractors LAWS Team members. The Baseline 
Backscatter profile is based upon data taken at JPL and 
WPL and represents the LAWS team’s best estimate of 
median values and variances. It must be noted that the 
profiles used to construct the baseline log normal 
distributions were the average of several 

100 individual profiles. Therefore, they do not contain 
the contributions of speckle to single shot backscatter. 
At this point in time, the contractors will need to 
consider the implications this has to the interpretation 
of LAWS performance. However, as soon as possible we 
intend to provide a single shot probabilistic 13 profile 
using the results of an analysis of RSRE data by Bowdle 
and Rothermel . 

2) Hosted an EosDIS Science Advisory Panel meeting 
addressing (among other topics) the needs to prototype 
the integration of LAWS with the baseline EosDIS. 

3) Participated in a GSFC workshop on using the NASA 
Climate Data System. Our specific interest is in 
using cloud climatologies based upon ISCCP and 
Nimbus-7 data sets to develop reasonable CFLOS 
statistics for LAWS as well as subvisual cirrus 
estimates based upon the average visual cloud properties. 

4) Received a funded extension (* $10K) to support 
the ongoing OSSE at FSU ( Krishnamurti ) . 


During the next period we will focus upon the following: 

1) providing a baseline velocity variance estimator 
for LAWS line-of-sight measurements. John Anderson, 
Bob Lee, Mike Hardesty and G.D. Emmitt are currently 
addressing this issue. 

2) modifying the Baseline Backscatter profile to include 
speckle statistics. 

3) providing a Baseline for av(z), dv/dz and correlated 
horizontal wind fields to complete the Baseline 
Atmosphere. The selection of <jv(z) is being done 
with consultation with other team members and active 
researchers in atmospheric turbulence. 
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4) preparing two papers for the Lake Tahoe meeting in 
February 1990 . 



BASELINE BACKSCATTER/ ATTENUATION PROFILE 


A probabilistic 0 profile, Figure 1, and tropical maritime 
attenuation profile, Figure 2, are provided to the contractors 
(GE and Lockheed) for use in their system models and for 
communicating their performances with selected LAWS 
configurations . 

The purpose of these profiles is to have a common basis for 
comparing the various trades that are being studied by not only 
the Phase A/B contractors but also the LAWS Science Team. 

The 0 profile is a smoothed version of a composite profile 
based upon ground-based observations at JPL and WPL. A consensus 
between several members of the LAWS Science Team (Bowdle, Post, 
Menzies and Emmitt) was reached regarding the following points: 

1) not enough data is available at 9.11 or 10.6 pm 
above 15 km to include this area in the baseline 
profile ; 

2) the general distribution of backscatter values 
around the median are log normal at all levels; 

3) both JPL and WPL data sets exclude the contribution 
of thin cirrus to the upper tropospheric backscatter; 

4) in situations without identifiable cirrus, the 
median 0 above 8 km is 3 x 10“ 11 m" 1 sr~ 1 (note: 
this value is not the lowest that could be argued 

to be consistent with the data collected at Mona Loa 
in 1988. Furthermore, this value is more repre- 
sentative of current conditions rather than those 
during the JPL/WPL monitoring period (since 84)) , 

5) the distribution of thin cirrus (t < 1.0) above 
8 km can be only estimated at this time. The 
value of 50% at levels above 14 km is a compromise 
between a value of ~ 30% seen from a ground 
perspective at JPL, Boulder and Hawaii and a value 
of 70-80% (tropics) that may be realized from a 
space perspective. The distribution of cirrus 
related backscatter is also assumed to be log normal; 
and 

6) performance implied by these profiles will be 
modulated by opaque clouds. For example, only 25-35% 
of all lidar shots in the tropics may reach, the 
earth’s surface due to clouds. 

The tropical attenuation profile is taken from LOWTRAN 7 
code and is provided to represent nearly 30-40% of the earth s 
surface (tropical maritime). It is the most severe mean profile 



and therefore should be used to "bracket" the performance. 
However, attenuation due to cirrus has not been included at this 

time . 


The software provided on disk is designed to return an 
answer to questions like the following: 

What % of the time will there be sufficient backscatter 

at 5 km to get a 5 dB SNR? 

In Figure 3, we show an example of profiles of probable 
performance using the baseline LAWS system as described on the 
attachment entitled "Baseline LAWS SNR Equation . 

Questions on the profiles and the software should be 
directed to Sid Wood or Dave Emmitt at SWA* 


Baseline LAWS SNR Equation 


% 


Several forms of the lidar SNR equation are available for 
use with the Baseline Atmosphere. After discussion with 
Hardesty, Bilbro, Menzies and others we have selected the 
following version: 


7t*c*ni *n 2 *n3 • n4 • J • D 2 *T*0*e" 2 / a(r)dr 

SNR = 

8 *hv • (R 2 + (.25-D • D/X) 2 ) 

c = speed of light (m/s) = 3.0 x 10® 

ni = heterodyne quantum efficiency = .40 

n2 = optical efficiency = .25 

n3 = beam shape factor = .46 

n4 = truncation factor = .70 

J = laser power (Joules) = 10 

D = mirror diameter (m) = 1.5 

t = pulse length (sec) = 6.67 x 10“ 6 

0 = backscatter ( m~ 1 sr _1 ) = input [ ~ 10 _1 °] 

e =2 way attenuation = computed [~ .1 - 1.0] 

hv = photon energy (J) = 1.88 x 10“ 20 

R = slant range (m) = computed 

= laser wavelength (m) = 9.11 x 10 -6 


X 
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TROPICAL MARITIME ATTENUAT I ON PROFILE: LOWTRAN 7 
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Figure 2. Average tropical maritime attenuation profile taken from LOWTRAN 7 code. While 

an average profile is shown here, it is recognized that there is a large variation 
about the mean in the real world. However, for trade studies this variation is no 
included. * 



BACKSCATTER PROBABILITY DIAGRAM FOR 5 DB SNR (9.2 uni) 
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PROBABILITY OF BACKSCATTER (*) GREATER THAN 
OR EQUAL TO THE BACKSCATTER NEEDED FOR 5 DB SNR 






Baseline LAWS LOS Velocity 
Error Estimates 


As with the lidar SNR equation, there are several radial or 
LOS velocity error estimates, or , that have been suggested for 
use with LAWS. While the Cramer-Rao Lower Bound may provide a 
limit to the extraction of a velocity estimate from a noisy 
signal, we have chosen the more conservative estimate based upon 
pulse pair autocorrelation processing of the Doppler signal. The 
following is derived from Eq. 6.22a in Doviak and Zrnic (1984). 

4 

V» ax = 
f 

t = 

w 

Vbw = 

Vatu = 

SNRw = 

Example : 

Given : 

SNRw = .079 (= -11 dB) 

W = .01 • ( .047 + 1 .0) • 5 = .01 
Or = .658 • (.112 + .200 + 160. 2)- 5 
or = 8.35 m s - 1 

See Figure A for a plot of SNRn vs o r . 


Or = 

X 


(2 x 1 ♦ 5 W + 16 x 2 W 2 /SNRw + 1/SNRw 2 ) • 5 
wavelength (m) = 9.11 x 10 -6 
maximum velocity measured = 50 ms -1 
sampling frequency = 2 • V«ax/X = 10.98 x 10 6 
pulse duration (sec) = 6.67 x 10 -6 

normalized frequency spread of return signal (m/s) 
1/Xf (Vbw 2 + V»t« 2 I 1 ' 2 

uncertainty due to pulse bandwidth (ms -1 ) = 

X 

2xt 

uncertainty due to turbulent eddies and wind 
shear within the pulse volume (m s -1 ) = 1.0 

W SNR 

SNRn = 5 dB = 3.162 
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A. 


We continue to focus our efforts on establishing baselines 
for evaluating various LAWS configurations. In particular, we 
are testing some LSM code that is used to generate simulated 
three-dimensional wind fields with prescribed shear, turbulence 
and mesoscale structures. We intend to present that simulation 
at the LAWS Team meeting in January. 

Considerable effort has been directed towards LAWS matters 
related to the Eos program’s science objectives and the EosDIS 
design activities. Many of the issues raised regarding expected 
LAWS standard data products have helped to define some areas for 
near term investigation. For example, the possibility of LAWS 
providing critical information on the presence and amount of low 
optical thickness cirrus needs to be carefully evaluated to 
ensure that appropriate attention is paid to calibration, shot 
management and signal processing options. 

Two papers have been submitted for presentation at the 
February meeting in Lake Tahoe on Optical Remote Sensing of the 
Atmosphere (see Appendix A). 

A memo was sent to the EosDIS Project Office advising them 
of a data masking plan for consideration by the LAWS team and the 
Eos project. A copy of the memo is attached along with a 
solicitation for comments. 

B. 


During our two-month no cost extension we will prepare a 
final report that will be submitted as a contractor’s report for 
NASA publication. 



November 8, 1989 


MEMO TO: 
MEMO FROM: 
SUBJECT: 


EosDIS Project Management 

G.D. Emmitt, EosDIS Science Advisory Panel Member 
Draft Ratipnale for LAWS Default Data Mask 


An underlying p 
prompt availability 
commonly used engine 
f ormats/ pro ject ions . 
after sensor transmi 
dedicated to having 
evaluated processing 
quality checking and 


rincipls of the EosDIS philosophy is the 
of remote sensing products expressed in 
ering units and universal 

Prompt is defined for EOS as a day or so 
ssion of the data. The EOS program is 
a full-up DIS with fully implemented and 
algorithms that will provide 1st order data 
sensor product generation* 


With the emphasis on prompt product generation, we can 
appreciate a second EosDIS principle - data Quality Assurance 
(QA) will, in effect, be carried out by the thousands of users 
each having different opportunities to expose data quality issues 
through different mixes of remotely sensed data, m situ data ana 
model outputs. The consequences of this QA by the user are o 
concern to the LAWS panel, particularly because of the 
exploratory nature of a space-based lidar facility. 


We recognize s 
rapid dissemination 
not all of the user 
to have the burden 
Worse yet are the c 
experienced in reco 
product as reliable 
resources in reconc 
inconsistencies . 


ome of the advantages and disadvantages of the 
of the EOS products. We also recognize that 
s of EosDIS products will appreciate or desire 
of QA placed upon them in the first instance, 
onsequences of users of EosDIS not adequately 
gnizing "bad" data - users accepting a NASA 
and therefore expending valuable time and 
iliation of data set conflicts and 


Another perspective on this issue is that of the scientis 
responsible for providing a standard data product. Regardless 
the accompanying documentation, error discussions and caveats, 


users will still associate a 
or group of scientists - i.e 
products* There is bound to 
scientists to contribute any 
producer . 


data set with a particular scientist 
, , they are perceived as "signature 
be a reluctance on the part of man> 
product not adequately QA ’ d by the 
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The LAWS team is considering the following as a more 
desirable option for the EosDIS product line - particularly in 
the first few years after launch. The concept is a Default Data 
Mask (DDM) . The DDM is simply a data quality mas p ace o 
the complete delivered data set so that a casual user or even a 
user that does not want to be involved in the QA of someone 
else ’ s data set will only receive those data that meet e 
highest quality standards. In other words, to get the entire 
s£lndard product, a user has to consciously remove the mask and 
accept the responsibility of using data not fully evaluated or 
endorsed by the product generating scientist. 


The advantages of the DDM option are. 


1 ) 


Standard products that have received^ 
possible within the "production mode 
delivered in their entirety to EosDIS 
"prompt availability" principle; 


the best quick QA 
time frame, are 
in the spirit of 


2) While having done their best, scientists (who know that gi 
time, better estimates of errors will be found f °f\ 
regions of marginal sensor performance) will be able 
to act more responsibly by masking those questionable 
data from the general user; 


3) The entire data set is available, however, to be used 
by those who are trained and prepared to recognize t e 
limits of sensor performance regardless of the plausi 
bility of the data; and 


4) The DDM can be slowly removed (or enlarged if necessary) 
as more is learned by the EosDIS community of users and 
confidence is established through validation and 
application . 


While this DDM approach may be most desirable for new 
technology" sensors like LAWS, other instruments with a space- 
based heritage (MODIS) may also find it desirable and an act 
scientific responsiblity to the broader spectrum of users being 
courted by the EOS program. 


It is requested that the EOS project provide the LAWS team 
with its reaction to this proposed way for a PI or Team Scientis 
to deliver a complete standard product to EosDIS while Providing 
some measure of default protection to the user who expects a QAd 
product . 
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Comments and suggestions should be directed to 

Chairman of LAWS Algorithm Committee 
G . D . Emm i 1 1 

Simpson Weather Associates, Inc. 

809 E. Jefferson Street 
Charlottesville, VA 22902 

and copied to: 

LAWS Team Leader 
Wayman Baker 
NOAA/NMC 

Development Division 
WWB Room 204 
Washington, D.C. 20233 



